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Abstract 

We explain the hnear polarization of the Ca ii infrared triplet observed close to the edge of the 
solar disk. In particular, we demonstrate that the physical origin of the enigmatic polarizations 
of the 866.2 nm and 854.2 nm lines lies in the existence of atomic polarization in their metastable 
^1-^3/2,5/2 lower levels, which produces differential absorption of polarization components (dichro- 
ism). To this end, we have solved the problem of the generation and transfer of polarized radiation 
by taking fully into account all the relevant optical pumping mechanisms in multilevel atomic 
models. We argue that 'zero-field' dichroism may be of great diagnostic value in astrophysics. 
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Spectropolarimetry provides key information on the physical conditions and geometry of 
astrophysical plasmas otherwise unattainable via conventional spectroscopy Q . The remote 
sensing of cosmical magnetic fields via the Zeeman effect is the most well-known example. 
In the presence of a magnetic field, the magnetic substates of atomic and molecular energy 
evels split, leading to a frequency shift of the (AM = ±1) and vr (AM = 0) transitions 
^. If the Zeeman splitting is a significant fraction of the spectral line width, the polarization 
of the different components does not cancel out and an observable polarization pattern can 
then be produced via the emission process, as well as through dichroism and anomalous 
dispersion phenomena. The observed line polarization then depends on the strength and 
orientation of the magnetic field vector along the line of sight. 

Yet there is a more fundamental mechanism producing linear polarization in spectral 
lines. In the outer layers of stellar atmospheres, where light escapes through the stellar 
surface, the atoms and molecules are illuminated by an anisotropic radiation field. The 
ensuing radiation pumping produces population imbalances among the magnetic sublevels 
of energy levels (that is, atomic polarization), in such a way that the population of substates 
with different values of |M| are different. This is termed atomic level alignment. As a 
result, the emission process can generate linear polarization in spectral lines without the 
need for a magnetic field. This is known as scattering line polarization 0]. However, 
light polarization components will also he differentially absorbed when the lower level of 
the transition is polarized 0, |^. Thus, the medium becomes dichroic simply because the 
light itself is escaping from it. Since this process has nothing to do with magnetic fields, 
we call it zero-field dichroism. This mechanism may be of great diagnostic value because in 
weakly magnetized astrophysical plasmas, like the 'quiet' solar atmosphere, the contribution 
of the transverse Zeeman effect to the linear polarization is normally negligible. In order for 
zero-field dichroism to be operative, the atomic polarization of the lower levels must survive 
depolarizing effects, such as those produced by elastic collisions. Moreover, weak magnetic 
fields modify the atomic level polarization via the Hanle effect 0], where the magnetic 
field strength B (in gauss) that produces a significant change is -B ~ 1.137 x 10^^/(tiife5'^) 
(tiife and being the level lifetime and Lande factor, respectively). Since the lifetimes 
of the lower levels of the transitions of interest are usually much larger than those of the 
upper levels, it is clear that diagnostic techniques based on this mechanism are sensitive to 
relatively low densities and magnetic field strengths. 
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FIG. 1: (a) Scattering geometry. Ray propagating along the direction Q, with polar angle 9. The 
reference polarization direction for Q > is along e^, i.e., parallel to the stellar surface, (b) 
Grotrian diagram of our atomic model for Ca ii indicating the total angular momentum of the 
levels and the spectral line wavelengths in nm. 

In this letter, we demonstrate that zero-field dichroism is the key mechanism underlying 
the linear polarization of the Ca ll infrared triplet detected in observations close to the edge 
of the solar disk, in very quiet regions far away from sunspots. The observed fractional 
polarization amplitudes of the Ca ll lines at 849.8, 854.2 and 866.2 nm are Q/J ~ 0.035%, 
0.13% and 0.12%, respectively |i2jj, the polarization signal in the 866.2 nm line being espe- 
cially intriguing. This is a ^03/2 ^Pi/2 transition whose upper level (^Pi/2) cannot be 
aligned. Therefore, emission following anisotropic radiative excitation cannot produce linear 
polarization in this spectral line and its detection has been considered 'enigmatic' 

Consider a plane-parallel stellar atmosphere without magnetic fields. For reasons of 
symmetry, choosing the reference polarization directions as in Fig. 1, the polarization state 
of a radiation beam can be described in terms of the intensity / and the Stokes parameter 
Q for linear polarization. The transfer of I and Q along a ray path element ds is described 
by the following equations: 



where e/ and eg are the corresponding emissivities, while rij and rjq are the absorption 
and dichroism coefficients, respectively. In a spectral line, the emissivities depend on the 
excitation state of the upper level of the transition. Conversely, rjj and rig are given by the 
excitation state of the transition's lower level. We have considered an additional contribution 




eg - VqI - ViQ , 



ei - r]il - r]QQ , 



(1) 



(2) 



to e/ and rji from a background continuum assumed to be unpolarized (i.e., 77^°"* and e™"* = 
rij°^^B,y, being the Planck function), which is a very good approximation towards the red 
part of the solar spectrum (see Isl]). For simplicity, we neglect stimulated emissions 15 1. 

We choose the quantization axis for angular momentum along the Z-axis of Fig. 1. Since 
the radiation field is axially symmetric around this direction, no quantum interferences 
are induced between magnetic sublevels and the excitation state of each atomic level is 
fully described by the individual populations (occupation probabilities) Nm of each of its 
magnetic sublevels. Equivalently, we can describe the excitation state of each level of total 
angular momentum, J, by means of the following components of the atomic density matrix 

Po{J)=y,{-^y~''^'^KTl\ \^\nm, K = 0,1,2,. ..,2J, (3) 

\M -M J 

which are just linear combinations of Nm- Because the radiation field in the assumed model 
atmosphere has no net circular polarization, we have an internal symmetry in our problem: 
Nm = N^M] therefore, only components pQ with K even can be different from zero. Of 
particular interest is Pq, which is (2J + 1)~^/^ times the overall population of the level, 
and the alignment coefficient p?,, which quantifies the population imbalance between the 

n n 

magnetic sublevels. The line emissivities can be expressed as [3|, HOI 

e'r = eop'o + eoWj^j,^i3p' - l)pl, (4) 
el'''' = eoWj^j,^il-p')pl, (5) 

where the pp ^^^1 Pq values are those of the upper level of the line transition under con- 
sideration, p = cos6' (see Fig. 1), and cq = (/iz//47r)yl„£0,^7V\/2 J„ + 1 (where z/ is the line 
frequency, A^i is the Einstein coefficient for the spontaneous emission process. A/" is the total 
number of atoms per unit volume, and (pi, is the absorption profile), wj^j^ is a numerical 
coefficient whose value depends on the total angular momentum of the upper and lower 
levels, taking values -2^/2/5, V2/10, and for the lines at 849.8, 854.2, and 866.2 nm, 
respectively. On the other hand, rji and rjq are given by identical expressions, but using the 
Pq and pg values of the lower level, r/o = {hv / An) B ^u't'x-N' \/2J + 1 instead of eo (with B^^ 
the Einstein coefficient for the absorption process), and Wj^j^ instead of Wj^j^,, where the 
values are —2^/2/5, V7/5 and V2/2 for the aforementioned lines, respectively. 
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It can be deduced from Eqs. (H}-© and the expressions of e/^g and r^/ g, that the emergent 
ractional hnear polarization at the core of a strong spectral line is approximately given by 
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where ctq = Po/Po fractional atomic alignment of the level under consideration. In this 

formula the ctq values are those at optical depth t = J rjids along the line of sight, where 
r ^ 1. Eq. (jH)) shows clearly that the observed fractional polarization in a given spectral 
line has in general two contributions, one from the atomic polarization of the upper level 
caused exclusively by emission from a polarized upper level (aQ^Ju)), and another one from 
the atomic polarization of the lower level ((7g(Ji)). As shown here, the latter one (caused by 
the selective absorption resulting from the population imbalances of the lower level) plays a 
key role in producing the linear polarization pattern of the Ca ll IR triplet. 

The atomic model in Fig. 1 is sufficiently realistic for modeling the Ca ll infrared triplet. 
Since calcium has no hyperfine structure, the only levels that may carry atomic alignment 
are the upper level ^P3/2 and the two metastable levels ^03/2,5/2- Therefore, nine density 
matrix elements are required to describe the atomic excitation: the overall population of 
each of the five levels (i.e., v^2XTTpq( J)), the atomic alignment of the two levels with 
J = 3/2, and the pi and pg values corresponding to the level with J = 5/2. These irreducible 
tensorial components of the atomic density matrix are governed by the statistical equilibrium 
equations, = dpQ{i)/dt = ~J2k' RiKK'Po' (i) + J2jK''^jiKK'Po' (j), where the relaxation 
matrix has the three contributions Rikk' = Rukk' + J2u RiuKK' + caused by decays 
toward lower levels, excitations towards upper levels, and depolarization by elastic collisions, 
respectively, while the transfer matrix has contributions caused by decays from upper levels 
{TuiKK') and excitations from lower levels {Tuxk') 13 • The relaxation rates read 

RiiKK' = SkK' + ^KK' Q^, (7) 
RiuKK' = ^KK' BiuJq + SkK' BiuJ^ + 5kK' Ciu, (8) 

and the transfer rates are 



TuiKK' — Y -^^J^^—^l^KK'PK'^ui + SkK' PK' Cui], (9) 

TuKK' = y -^-j-^^l^KK' Pk' BuJq + qxK' BuJq + 6kk' Pk' Ca], (10) 
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TABLE I: Numerical coefficients of Eqs. (j7j)- (|lU() calculated for the 5-level Ca ii model of Fig. 1. 
Pq = 1 for all transitions. For the 396.9 nm line (J/ = Ju = 1/2) any other coefficients are zero. 
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where p, q, and s are numerical coefficients that depend only on the quantum numbers of 
the levels involved in the transition (see Table 1). Since the radiation field is axisymmetric 
the absorption rates depend only on the line-integrated mean intensity (Jq), and on 



Jo = f ^ d/^'0.[(3/x'2 - 1)/ + 3(1 - 



It is important to note that in stellar atmospheres the polarization degree is small (Q/J ^ 1), 
and the main contribution to Jq results from the intensity average. We call A= Jq/Jq the 
'degree of anisotropy' of the radiation field. It determines the fractional atomic polarization 
that can be induced by optical pumping processes. In the solar atmosphere ^ ^ 1, hence 
ctq ^ 1. CoUisional rates Cij due to collisions with an isotropic Maxwellian distribution of 
electrons are modeled in analogy to the excitation by an isotropic unpolarized radiation field. 
We have used the inelastic coUisional rates for Ca ll tabulated in [ll| , while the depolarizing 
rates D^^^ due to collisions with neutral hydrogen atoms have been estimated according to 
Ref. Q- 

Figure El shows results of the self-consistent solution of the radiative transfer and statistical 
equilibrium equations in a semi-empirical model of the solar atmosphere. The numerical 
solution has been obtained applying the iterative methods outlined in Ref. P]. The figure 
indicates how the degree of anisotropy corresponding to the K-line and to the 854.2 nm line 
vary with the K-line optical depth. These spectral lines are representative of the radiation 
fields in the UV doublet and IR triplet, respectively. The figure shows also the behavior 
of the fractional atomic alignments ctq of the 2nd, 3rd and 5th levels. Note that they have 
sizable values, even at the depths in the stellar atmosphere model where the line optical 
depths are unity along the line of sight. 

The atomic alignment of the 5th level is governed by the radiation field in the K-line, with 




FIG. 2: Variation of Pg/Po "^^^ ™^ ^^'^ levels of our Ca ii model as a function of the K-line 

optical depth calculated along the vertical direction (solid and dotted lines). Dashed lines show 
Jg/Jg in the K and 854.2 nm lines. Arrows indicate the layer in the FAL-C model atmosphere 
where, for each of the three IR lines, r = 1 along the p = 0.1 line of sight. 
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FIG. 3: Emergent fractional polarization of the Ca ii infrared triplet at ;U = 0.1 in the FAL-C semi- 
empirical model of the solar atmosphere (solid lines). The shapes of the calculated Q/I profiles and 
their relative line-core amplitudes are in good agreement with the reported observations (see [3|). 
Dotted lines: Q/I profiles if the /orif/- /iued metastable '^^3/2,5/2 levels were completely unpolarized. 



optical pumping in the lines at 849.8 nm and 854.2 nm playing only a marginal role. In fact, 
the height variation of cTq(5) closely follows that of the degree of anisotropy in the K-line 
throughout the whole atmosphere, with crg(5) ~ V2/2 {Jq/ Jq)(K). On the other hand, there 
are two optical pumping processes able to introduce atomic alignment in the lower levels 
of the Ca ll IR triplet (the metastable levels): repopulation pumping, which results from 
the spontaneous decay of the polarized 5th level; depopulation pumping, which occurs when 
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some lower state sublevels absorb anisotropic light more strongly than others. Since the 
infrared triplet lines are much weaker than the ultraviolet H- and K-lines (^7^|™ey'^cont ^ j^gs 
for the former, while 7^i.™<5^^cont ^ j^qT f^j, ^^le latter lines), the IR line photons escape from 
deeper atmospheric regions, where the Jq / Jq values of the IR line transitions are much larger 
than the (negligible) cTq value of the 5th level at such depths (see Fig. Consequently, 
at such atmospheric depths the long-lived metastable levels 2 and 3 are mainly aligned 
by depopulation pumping rather than by spontaneous emission from the upper level with 
J = 3/2. 

The solid lines in Fig. IHl show the emergent fractional polarization of the Ca ll infrared 
triplet calculated at /i = 0.1. Dotted lines show the emergent fractional polarizations if 
the metastable levels were unpolarized and dichroism negligible. Since the upper level of 
the 866.2 nm line cannot be aligned, the Q/I signal in this transition is exclusively caused 
by dichroism [see Eq. ©]. The linear polarization observed in the 854.2 nm line can, in 
principle, be caused by both, the emission events from the polarized upper level and to the 
differential absorption by the polarized lower level. However, |(To(^D5/2)| > |crQ(^P3/2)| at 
r(854.2) ^ 1 (see Fig. E)), while > wj^j^ in Eq. (jH}. Therefore, the observed Q/I 

signal is de facto also mainly caused by dichroism in the solar atmosphere. The 849.8 nm 
line forms significantly deeper in the solar atmosphere, where the fractional alignment of its 
upper and lower levels are close to zero, hence its much lower polarization signal, whose sign 
turns out to be very sensitive to the assumed atmospheric model. 

Note that the calculated Q/I amplitudes of the 854.2 nm and 866.2 nm lines (see Fig. 
are roughly twice the values reported in This can be due either to a simplified atmo- 
spheric modeling (the highly inhomogeneous solar atmosphere may not be well represented 
by a one-dimensional semi-empirical model), and/or to the existence of depolarizing mech- 
anisms. Either case emphasizes the great diagnostic potential of the scattering polarization 
observed in the Ca ll IR triplet. For example, coUisional depolarizing rates five times larger 
than the nominal values used here only slightly affect the ensuing Q/I amplitudes. On the 
other hand, in the presence of a microturbulent magnetic field with a strength as low as 0.01 
gauss, the atomic polarization of the metastable ^03/2,5/2 levels are significantly reduced 
and hence the emergent linear polarization in the 854.2 nm and 866.2 nm lines decreases to 
values similar to those observed (Hanle effect). 

Zero-field dichroism may also be operating in other astrophysical objects (e.g., accreting 
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systems) and should be fully taken into account when interpreting spectropolarimetric ob- 
servations in other spectral lines besides the Ca ll IR triplet itself, whose polarization has 
been observed recently in supernovae 
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For a blackbody at the effective temperature of the Sun (T ~ 5800 K), the number of photons 
per mode even at 900 nm is n ~ 6 x 10^^ <ti 1, and stimulated emission can be safely neglected. 
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